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ABSTRACT 

The dynamics of two initially unmagnetized relativistic counter-streaming homogeneous ion-electron 
plasma beams are simulated in two dimensions using the particle-in-cell (PIC) method. It is shown 
that current filaments, which form due to the Weibel instability, develop a large scale longitudinal 
electric field in the direction opposite to the current carried by the filaments as predicted by theory. 
East moving ions in the current filaments decelerate due to this longitudinal electric field. The same 
longitudinal electric field, which is partially inductive and partially electrostatic, is identified as the 
main source of acceleration of electrons in the current filaments. The transverse electric field, though 
larger than the longitudinal one, is shown to play a smaller role in heating electrons, contrary to 
previous claims. It is found that, in ID, the electrons become strongly magnetized and are not 
accelerated beyond their initial kinetic energy. Rather, the heating of the electrons is enhanced by 
the bending and break-up of the filaments, which releases electrons that would otherwise be trapped 
within a single filament and hence slow the development of the Weibel instability (i.e. the magnetic 
field growth) via induction as per Lenz’s law. In 2D simulations electrons are heated to about one 
quarter of the initial kinetic energy of ions. The magnetic energy at maximum is about 4 percent, 
decaying to less than 1 percent by the end of the simulation. Most of the heating of electrons takes 
place while the longitudinal electric field is still growing while only a small portion of the heating is 
a result of subsequent magnetic field decay. The ions are found to gradually decelerate until the end 
of the simulation by which time they retain a residual anisotropy less than 10 percent. 


1. INTRODUCTION 

Collisionless shocks forming in astrophysical environ¬ 
ments are believed to be mediated by electromagnetic in- 


sta lpilities dKennel and Sagdeev|1967|j Parker|1961| |Eich- 
I^|1979 Blandford and Eichlerj 19871 ) in which ions are 


scattered by the resulting magnetic field fluctuat ions. In 


particlular, the Weibel instability (Weibel 1959), which 
causes fast growth of strong magnetic held at small-scale 
in anisotropic plasma flow, has received much attention 
as the main isotropization mechanism that leads to shock 
transition in free-streaming ejecta from violent astro- 


physical events ( 

Medvedev and Loeb|1999 Wiersma and 

Achterberg||200^ 

l[ Lyubarsky and Eichler 2UU6[ |Achter- 


berg and Wiersma 120071 

IBretl 

120091 Yalinewich and 

Gedalin 2010| Shaisultanov 

■ et a 

/. |2012). As the domi- 


ion gyroradius, which renders them inefficient scatterers 
of ions, the longstanding question about its role in col¬ 
lisions shocks ha s been how well it competes with other 
mechanisms (e.g. Galeev et a/.|1964l Blandford and Eich¬ 
ler 1987 Lyubarsky and Eichler ' ^06 [ ). Clearly mecha¬ 


nisms that require a preexisting magnetic field are ques¬ 
tionable when the magnetic field is weak. On the other 
hand, a slow shock mechanism may actually suppress 
a faster mechanism because it creates a broader shock 
transition, and thus has greater ’’reach” upstream of the 
shock. 

The role of Weibel instability in forming the shock 
transition in weakly magnetized plasmas has recently 
been established through several numerical experiments 


(iNishikawa a/.|2QQ3 |Frederiksen^^^^.|2QQ4[ 

et a/.|2UU5t bpitkovsky 2UU8t|Keshet et a/.|2UU9 


JN umeri- 

cal simulations of relativistic shocks have shown that the 


pUU! 


Nishikawa 


relatively less energetic upstream electrons are signifi¬ 
cantly heated to energy comparable to the energy of ions 
as they cross the foreshock, which appears to be in agree¬ 
ment with the electromagnetic observations of supernova 
remnants and gamma-ray bursts afterglows, where the 
high energy radiation is believed to be synchrotron radi¬ 
ation originating from gyration of high energy electrons 



in the foreshock region, both in nature and in the numer¬ 
ical simulations of collisionless shocks, and whether the 
Weibel instability-induced magnetic field should persist 
over astrophysically significant length scales behind the 
shock. 

Numerical simulations suggest presence of large scale 
electric field in and around the current filaments forming 
just ahead of the shock and ions are found to decelerate 
due to this electric field. The same electric field that de- 


celerates ions should also accelerate electrons dBlandford 

and Eichler 

19871 

ILyubarsky and Eichler| 

20061 Ifiedalm 

et al. 2008 

2012| 

). However, the detail 

is of the elec- 


tric held and the acceleration mechanism have not been 
worked out. 


In order to understand the heating of electrons in fore¬ 
shock region of relativistic shocks we simulate, using the 
kinetic PIG method, development of the Weibel instabil¬ 
ity in two relativistic counter-streaming plasma beams, 
which resembles the precursor of a relativistic shock, but 
simpler and more idealized and therefore more suitable 
for resolving certain fundamental questions. The beams 
are taken to be homogeneous and fully interpenetrating 
at the beginning, so the time development of the instabil- 
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ity in our simulation imitates the spatial development of 
the instability in the foreshock regions, where later times 
in our simulation correspond the regions closer to the 
shock fronts. Our quantitative analysis suggests that the 
most of the heating of upstream electrons is due to the 
longitudinal electric field. The transverse electric field, 
though much stronger than the longitudinal, has negligi¬ 
ble effect on net acceleration of electrons. 

We have also inserted virtual test particles into the 
code to act as numerical probes. They assist in the di¬ 
agnostics of the mechanisms that are at work, and this 
will be discussed below. 

2 . NUMERICAL SIMULATION 


We use 

parallel version of the PIC code TRISTAN 

(Buneman 

|1993 Spitkovsky 

2005) to simulate relativis- 


tic counter-streaming beams of ions and electrons in two 
dimensions. The simulation is initialized by placing ion- 
electron pairs at uniformly chosen random locations in 
a two dimensional square box in x-y plane. Initially ion 
and electron in each pair are moving in opposite direc¬ 
tions along the x-axis with Lorentz factor 70 . Half of the 
total number of pairs have ions moving along the posi¬ 
tive x-axis while an equal number of another half pairs 
have ions moving along the negative x-direction, hence 
creating two oppositely streaming neutral plasma beams 
of equal intensity, where each beam has a net current. 
The initial condition ensures that the simulation is ini¬ 
tially charge and current free and that Maxwell equations 
are satisfied over scales where there are equal numbers of 
forward and backward moving electron-ion pairs, but at 
the grid scale, local fluctuations in the number of forward 
and of backward movers imply a net small scale current. 
This noise level, however, is reddened by smoothening 
the grid scale current. As the electrons and ions in the 
pairs separate from each other, the electric and mag¬ 
netic field grow from the noise in the current generated 
by streaming electrons and ions. Since we simulate an 
initially unmagnetized plasma, only out-of-plane mag¬ 
netic field Bz and in-plane electric fields and Ey are 
excited as the Weibel instability set in. We impose pe¬ 
riodic boundary condition in both x and y directions for 
both particles and fields. 

Here we discuss results mainly from the two largest 
simulations we have attempted in terms of physical size 
and evolution time. Ion to electron mass ratio rrii : rrie 
for the two reported simulations are 16:1 and 64:1, and 
are henceforth referred to as Miq and M 64 , respectively. 
All the figures in this paper are for ion to electron 
mass ratio of 64:1, i.e. simulation M 64 , unless otherwise 
stated. The physical sizes of the box for Miq and Mq/^ 
are 500 c/cjpexSOO c/uOpe and 1000 c/cjpexlOOO cjujpe^ re¬ 
spectively, where c is the speed of light in vacuum and 
ujpe = Y^47rnoe^/me7o is the initial electron plasma fre¬ 
quency, where e is the charge of an electron and no is the 
initial number density of electrons, or by charge neutral¬ 
ity, of ions. Both simulations were resolved to 1/lOth of 
the initial electron skin-depths. Initially there are 32(8) 
particles per unit cell of the simulation box for Mi 6 (M 64 ) 
and the simulation was evolved for 500 (1000) plasma 
time uJpQ. Initially the ions and electrons are moving 
along the x-axis with initial Lorentz factor 70 = 10 in 
the both cases. 


3. FILAMENTATION AND ELECTRON HEATING 

At the very beginning of the simulation the counter¬ 
streaming electrons, which are relatively lightweight, as 
compared to ions, and can relatively easily be deflected 
by the magnetic perturbations, are subject to the Weibel 
instability. Current due to the streaming electrons gener¬ 
ate magnetic field and the electrons moving in the same 
direction are herded into the filaments by self created 
magnetic field. Growing current in the filaments induces 
an electric field opposite to the direction of current in 
accordance with Lenz’s law. The induced electric field 
slows down the streaming electrons in the filaments and 
are scattered by the strong magnetic field around the fila¬ 
ments. During this stage relatively heavier ions continue 
streaming almost unaffected. As the electron Weibel in¬ 
stability stage ends in about 20 uj ~^, electrons are nearly 
thermalized to about their initial streaming kinetic en¬ 
ergy noTYie^jo — l)c^, creating a nearly isotropic electron 
background for still streaming ions. 

As the electron filamentation stage ends, relativisti- 
cally counter-streaming ions undergo the Weibel instabil¬ 
ity and likewise current carrying ions are separated into 
long current carrying filaments along the initial stream¬ 
ing direction. The current filaments are non-stationary 
structures with slight bending in the transverse direction 
(figure [B , which move along with the current carrying 
ions with the strength of magnetic field being small near 
the filaments and large between two adjacent filaments. 
Again, the inductive electric field develops in the fila¬ 
ments opposite to the current due to the streaming ions. 
In addition to the inductive electric field along the fila¬ 
ments an electrostatic electric field develops around the 
filaments in the direction transverse to the streaming di¬ 
rection due to excess of positively charged ions in the 
filaments. The scale of this transverse electrostatic field 
is limited by the electrons that cloud around the fila¬ 
ments and remains ^ cjujpe^ which grows with time as 
the electrons are accelerated to higher energies. In Fig¬ 
ure we show the details of the local magnetic field and 
the local current due to ions and electrons in a small spa¬ 
tial part of the simulation box. The apparent correlation 
between the motion of ions and electrons (electrons move 
opposite to their current) suggests that the electrons are 
accelerated in the direction of streaming ions. 

The spatial distribution of electrons closely follows the 
distribution of ions. Charge neutralizing electrons clump 
in rather positively charged filaments and stream along 
with the ions. The electrons in the filaments move along 
the inductive electric field created by the streaming ions 
and consequently get energized. Ions on the other hand 
move against the electric field induced by themselves and 
hence lose their energy. Since the current filaments are 
non-stationary and turbulent at very small scale, and 
so are the field structures, the local details of energy 
gain and loss can deviate significantly from the one men¬ 
tioned above and hence an statistical approach is needed 
to understand the systematic acceleration of electrons. 
We therefore compute spatial average of various quanti¬ 
ties and correlates to quantify the presence of systemic 
heating of electrons at scales much larger than their skin 
depth. 

During the linear stage of the instability not all the 
energy that is lost by ions goes into heating electrons. 
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Figure 1. Left panel: Structure of magnetic fields at 225 ojpe ■ 
The magnetic field is still growing due to the Weibel instability 
and forms filamentary structures parallel to the streaming direction 
(referred to as filamentary phase). Most of the heating of electrons 
takes place during the filamentary phase. Right panel: Structure of 
magnetic field at 650 ujpe ■ Transverse size as well as the bending of 
filaments grow with time and eventually filaments get disoriented. 
Here, and henceforth, Bz and all other electromagnetic fields (i.e.. 
Ex and Ey) are normalized to (47mo(7o —unless 
otherwise specified. 



Figure 2. Local magnetic field, and the current due to electron 
and ions are shown in a spatial patch from the simulation Mq/^ ( 
rrii : me=64:1) at time 65 ujpe • Direction of black and red arrows 
indicate direction of the local current due to the ions and electrons, 
respectively, at the tail of the arrows and length of the arrow is in 
proportion to the magnitude of the corresponding vector quantity. 
The local currents shown here are computed by taking the aver¬ 
age velocity of all particles within a box of size ^Acuope x ^Acuope 
centered at the tail of arrows. Inductive electric field opposes the 
current due to fast moving ions in the current filaments and is 
responsible for deceleration of ions as well as acceleration of elec¬ 
trons. 


A significant part of the total kinetic energy goes into 
the electromagnetic fields which are responsible for me¬ 
diating collective interactions between the charged par¬ 
ticles. In figure we show how much of initial kinetic 
energy goes into different components of the electromag¬ 
netic fields. During the ion filamentation stage, which is 
mainly magnetic in nature, energy in the magnetic fields 
reaches few percent of the total initial kinetic energy. 
Longitudinal electric field Ex is the weakest component, 
but as shown later, is the most important for the net 
transfer of kinetic energy from ion to electrons. The lon¬ 
gitudinal electric field Ex has a curl in the direction of 
magnetic field which predominantly contributes to the 
growth of the magnetic field. In other words, the induc¬ 
tive electric field Ex accompanies the growth of the mag¬ 
netic field and accelerates electrons hence leaving it the 
weakest amongst all growing electromagnetic field com¬ 
ponents. As the non-linearity sets in and the filaments 
are bent and broken. Ex and Ey both converge to a sim¬ 
ilar value since disruption and disorientation of current 



Figure 3. Temporal variation of the mean square (averaged over 
the physical domain of the simulation) of the electromagnetic field 
components Ex (dashed red), Ey (dotted dashed green), and Bz 
(solid blue) are normalized to 47mo(mi +me)(7o — l)c^. Amongst 
the electromagnetic fields, the magnetic field gains most of the ki¬ 
netic energy lost by the ions to the fields, which holds positively 
charges ions in the filaments together against the electrostatic elec¬ 
tric field, followed by mostly electrostatic transverse electric field 
Ey and longitudinal electric field Ex- 


filaments wipe out any preference in the direction of the 
electric field in the x-y plane (figure [^. 

3.1. electron heating estimate 

The inductive longitudinal electric field in the direction 
opposite to the current in plasma features in the linear 
theory of the Weibel instability itself. In fact it is the 
curl of this longitudinal electric field which accompanies 
the growth of magnetic field in accordance with the Fara¬ 
day’s law of induction. Presence of a large scale electric 
field in a conducting plasma naturally predicts transfer of 
kinetic energy from ion to less energetic electrons. Here 
we estimate the net energy gained by electrons due to 
the inductive electric field in a one dimensional linear 
theory of the Weibel instability. In the two dimensional 
case simulated here electromagnetic field fluctuations due 
to the electrostatic wave modes parallel to the x-axis 
(streaming direction) and the oblique modes also become 
comparable to the transverse modes, and eventually lead 
to the disruption of current filaments (Shaisultanov et al 
2012). Here we consider the heating solely due to the 


transverse modes with the wave vector k along the y- 
axis and electric field E along the x-axis. We assume 
homogeneity along the x-axis, that is averaging out the 
fluctuations along the x-axis, hence reducing the problem 
to one spatial dimension. 

Let Px{y) be the x-momentum of an electron at any 
given location y. The electric field Ex{y) at any given 
y determines the rate of change of x-momentum of an 
individual electron located at y, i.e. dpx{y)/dt = eEx{y). 
Averaging the x-momentum of all electrons at a certain 
y gives. 


dp:c{y)/dt = eE^{y) (1) 

where p^^iy) = J Pxfe{y,P,t)d^p/ J fe{y,P,t)<rp with fe 
being the distribution function of electrons which is as¬ 
sumed to depend on y only. The electromagnetic fields 
Ex and Bz , as well as the mean x-momentum of electrons 
px , can all be written as sum of sinusoidal modes of vari¬ 
ous wavelengths, but with Bz quarter a wavelength phase 
shifted from Ex and Px^ since Bz is maximum between 
two adjacent current filaments but Ex and Px peak in the 
current filaments. We first consider a single sinusoidal 
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mode of wave number k and write Bz = B^^hi{ky)^ 
where B^ is time dependent amplitude of the magnetic 
field Bz in the mode of wave number k. From Maxwell’s 
equations we relate the electric and magnetic field as 
kcEx = cos{ky)dB^/dt. Decomposing the linear equa¬ 
tion!^ into Fourier components and then substituting for 
the electric field we find 

dp^/dt = {ejkcPjdB^ jdt (2) 

where p\ is the amplitude of average x-momentum px 
in the mode of wave number k. Equation enables 
us to express instantaneous amplitude of average three- 
momentum fluctuation in terms of the amplitude of mag¬ 
netic fluctuation as 

pl = eB’^/kc + pM (3) 

Equation [^is of central importance: It implies that the 
gyroradius m an electron, in the limit that Px » Pxo^ 
is roughly the wavelength times l/27r. As we shall see 
below, this is confirmed by the simulations. The impli¬ 
cation is that the field growth is not stopped by the mag¬ 
netization of the electrons, as conjectured by Lyubarsky 
and Eichler ( 2006), because the electrons are at all times 
only marginally magnetized. That the energy of the elec¬ 
trons keeps pace with the field growth, so that they never 
get highly magnetized, is the essential reason the elec¬ 
trons get heated as much as they do. 

It is also shown that the marginal magnetization on the 
other hand, is enough to suppress the electron heating 
in one dimension, because the filaments remain exactly 
straight, and the electrons are accelerated exactly along 
the filaments. By contrast, in more than one dimension, 
the existence of oblique modes allows the filaments to 
bend, and it becomes harder for an electron to remain 
within a single filament. The countercurrents of the elec¬ 
trons are then less likely to cancel those of the ions and 
the field growth can proceed to larger length scales. 

At any instant, the rate of change of average kinetic 
energy of an electron Ue = ^e(7e — l)c^ due to the lon¬ 
gitudinal electric field Ex is given by the spatial average 
(y-average) of eExVx- For a single mode, defining as 

g -average of cp^ we obtain (after squaring the equation 
and then taking the derivative with respect to time) 

U^dU^ldt ^ \e^c^l{kcf\B^dB^ldt (4) 

The linear theory of Weibel instability predicts exponen¬ 
tial growth of the amplitude of magnetic field fluctuation 
B^ at all length scales. However, the maximum strength 
of magnetic field at any given length scale is constrained 
by the available current in the plasma at that length 
scale. Consequently, the growth of magnetic field am¬ 
plitude B^ at a length scale 1/k reaches a saturation 
amplitude which is inversely proportional to the 

wave number *0 Since the amplitude of magnetic field 
B^ at smaller scales saturates first, while the magnetic 
field at larger scales are still growing, the net magnetic 
field at any instant is dominated by the magnetic field in 


^ the maximum achievable net current in the filaments of trans¬ 
verse size 1/k is limited by the density of streaming ions in that 


length scale, which is ~ noec/A:, implying that oc 1/A: (Kato 

|2^|G^lin et ol|20lo| 




Figure 4. Left Panel: average (amplitude of magnetic field 
in mode of wavenumber k) as a function of wave number k at time 
50, 150, and 250 oj^e are shown by red, blue, and green curves, 
respectively. The spectrum is obtained by first computing one 
dimensional discrete Fourier transform of magnetic field in several 
slices of the simulations box along the transverse direction and then 
taking average of the magnetic field amplitudes in any given mode 
k over all these slices. Right panel: evolution of average amplitude 
of transverse magnetic fluctuation for 27r/k = 10, 2,1, O.bcujpp 
are shown by red, green, blue, and black curves, respectively. The 
magnetic field in any given mode grows exponentially and then 
slowly decays. 




Figure 5. Transverse spectrum, as in the left panel of figure^ of 
longitudinal electric field Ex and transverse electric field Ey as a 
function of wave number k at time 50, 150, and 250 ouPe are shown 
by red, blue, and green curves, respectively. 


the largest transverse scale achieved during the Weibel 
instability by that time. In Eigurefflwe show the time 
evolution of transverse spectrum of the magnetic field 
structure in the simulation which confirms the exponen¬ 
tial growth and saturation of the magnetic field at length 
scales smaller than a few times the skin depth of ions. In 
figure!^ we show the transverse spectrum of Ex and Ey 
as welT It shows the exponential growth and satura¬ 
tion of the transverse electric field, and that the trans¬ 
verse electric field peaks at half the wavelength magnetic 
field peaks, since the transverse electric field is symmet¬ 
ric with respect to the current filaments. The transverse 
spectrum of longitudinal electric field in figure does not 
show saturation because of growth of electrostatic par¬ 
allel modes, as shown later, which contributes to further 
growth of the longitudinal electric field . 

The heating due to the longitudinal electric field Ex is 
significant only until the filaments prevail. As the trans¬ 
verse scale of the filaments reach ions skin-depth, bend¬ 
ing of the filaments in the transverse direction becomes 
significant. The current filaments start to break, get ran¬ 
domly oriented in x-y plane, and the rate of heating of 
electrons is substantially reduced. In order to obtain an 
estimate for the total energy acquired by the electrons by 
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the end of the filament at ion it suffices to integrate equa¬ 
tion]^ with respect to time for the smallest wavenumber 
kmin achieved during the Weibel instability with the elec¬ 
tron heating in progress, and integrating it up to the time 
saturation is achieved in this mode, that is to say until 
the heating due to this mode is significant. Following 
these prescriptions for the net heating estimate we get, 

« TOi7oc^(wpi/fcmmc)^(Bp™^^/87rno) (5) 

where t// is the average energy of an electron towards 
the end of the heating process during the linear phase. 

As compared to the gyro-frequency of electrons, the 
transverse scale of filaments grow relatively slowly with 
time. Electrons gain energy in the current filaments and 
their Larmor radii get enlarged. However, at any instant 
during the linear stage of the Weibel instability the en¬ 
larged Larmor radii of electrons due to inductive electric 
field in the current filaments cannot substantially exceed 
the transverse size of the current filaments. This limita¬ 
tion on electron Larmor radii is due to the fact that for 
an electron that has Larmor radius much larger than the 
spacing between two adjacent filaments, energy gained 
in one filament is lost in the neighboring filaments which 
carries opposite current and hence directed in op¬ 
posite sense. Additionally, as suggested by equation]^ 
the energy imparted by the inductive longitudinal elec¬ 
tric field in the largest length scale is just enough to 
keep the Larmor radii of electrons at any instant about 
l/27r times the transverse size of the current filaments. 
As the electrons get accelerated to higher energies their 
skin-depth also increases in the same pro^rtion. In¬ 
deed, as observed in the simulations (figure]^, the ratio 
of the Larmor radius to the instantaneous skin depth 
of the electrons, which is also about the l/27r times the 
transverse size of the current filaments, is nearly constant 
during the linear stage of the Weibel instability. It sug¬ 
gests that during the filamentation stage the inductive 
heating of electrons proceed such that nolle ^ B‘^IS7t 
remains satisfied. 

The growth of magnetic field and transverse size of the 
filaments continue until the transverse scale of the fila¬ 
ments becomes comparable to the ions skin depth c/ujpi 
after which the filaments are disoriented and the heating 
of electrons due to longitudinal inductive field is sub¬ 
stantially reduced (in the two dimensional case electrons 
continue to gain energy due to decay of magnetic field, 
but at a much reduced rate). It suggests that for the 
purpose of estimating the total heating due to inductive 
electric field we can take ckmin ^ This along with 
the condition that during the heating nUe ^ sug¬ 

gest that I// < rriijoc^ (from equation]^, that is to say 
the electrons are accelerated to energy comparable to the 
energy of ions during the linear stage of the Weibel insta¬ 
bility due to the inductive longitudinal electric field. In 
the following we present results from the PIC simulation 
which confirms that the heating of electrons due to the 
inductive longitudinal field is indeed significant, though 
there is additional heating due to the electrostatic modes 
which appear in the two dimensional analysis and is es¬ 
sential to the heating of electrons. 

3.2. Net work done by the eleetrie field 




pe pe 

Figure 6. Left Panel : The mean Larmor radii of electrons as well 
as test particles (charged particles that do not participate in the 
plasma dynamics but respond to the local electromagnetic field) 
with charge to mass ratio 0.1 than that of the electrons are shown 
in units of instantaneous skin depth of electrons. The red curve 
shows Larmor radius to electron skin depth ratio for electrons, 
while black, blue, and green curves show the same for test parti¬ 
cles with initial kinetic energy 0.01, 0.1, and 10 times that of the 
electrons. During the linear stage of the instability Larmor radius 
of electrons as well as the test particles gaining energy from the lon¬ 
gitudinal electric field remains about the instantaneous skin depth 
of the electrons. Right Panel: Temporal variation of mean kinetic 
energy of the electrons and test particles in units of yorUeC^. Color 
code for the curves is same as in the left panel. Energetic parti¬ 
cles that have Lamor radii much larger than the spacing between 
the filaments, such as the test particles with initial kinetic energy 
(70 — l)meC^ (shown in blue) do not get efficiently energized by 
the longitudinal electric field. 
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Figure 7. Evolution of mean kinetic energy of electrons Ke 
(dashed red ) and ions Ki (solid green) are shown for two dif¬ 
ferent ion to electron mass ratios (normalized to the initial kinetic 
energy electrons K^o) . Electrons heating due to the Lenz electric 
field is significant until the ion current filaments are disrupted. By 
the end of the simulation electrons are heated to about 1/3 of the 
ions kinetic energy in both cases consider here. 


While the magnetic field scatter the charged particles 
by altering their trajectories, the only source of acceler¬ 
ation of the plasma particles is electric field in the lab 
frame. From our PIC simulations we separate out work 
done by the two orthogonal components of the electric 
field, namely and Ey^ on electrons as well as on pro¬ 
tons. As evident from Figure]^ in both the simulations 
reported here the change in total kinetic energy of elec¬ 
trons and protons is apparent to be mainly due to the 
longitudinal electric field E^. The net work done by the 
transverse electric field which is rather stronger than 
the longitudinal electric field (Figure]^, is vanishingly 
small El 

^ The implication is that the net energization of electrons due 
to electro ns falling into the grow ing a nd merging filamen ts, as sug¬ 
gested by |Hededal et al] (|2004| and |Spitkovsky| (|2008|) , is rather 
small during the hlamentation stage ot the Weibel instability. 
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Figure 8. The rate of energy gained (lost) by electrons (protons), 
as a function of time, due to Ex and Ey, i.e. ExVx and EyVy, re¬ 
spectively (in units of c^^dTryorUe), are shown by solid and dashed 
red (black) curves, respectively . The quantities ExVx and EyVy 
are the mean values of ExVx and EyVy, respectively, computed for 
the sample of particles (few percent of the total number of particles) 
which were initially homogeneously distributed over the simulation 
box. The left and right panels are for ion to electron mass ratios 
16:1 and 64:1, respectively. Though the transverse electric field Ey 
is much larger in magnitude than the longitudinal electric field Ex 
( figure]^, its contribution to the electron acceleration and ion de¬ 
celeration (dashed lines) is negligible as comparable to the heating 
due to Ex- 


In figure we show the time evolution of the total ki¬ 
netic energy of electrons and ions. It shows that most 
of the energy exchange between ions and electrons takes 
place during the filamentation stage of the instability, 
when magnetic field is still growing (figurein strength. 
We have run the simulation for long enough to capture 
most of the heating phase and long after the rapid heat¬ 
ing phase ends, electrons are found to have acquired a 
substantial part of the ions’s kinetic energy, in agree¬ 
ment with the estimate presented above. Even long after 
the breaking of filaments, ions are not completely ther- 
malized and continue to lose energy and the heating of 
electrons continues, though at a much lower rate. 

Bending of the current filaments is apparent from the 
very beginning of the Weibel instability. The bending 
of current filaments results in mixing of the longitudi¬ 
nal and transverse components of the electric field. Two 
dimensional linear theory of the Weibel instability pre¬ 
dicts the growth of waves with wave vector at oblique 


angle with respec t to the streaming direction (Shaisul- 
tanov et al. |2012|), and we suggest that this can account 
for the bending of the filaments. The filaments are also 
susceptible to the Buneman instability, which leads to 
the growth of resonant waves parallel to the streaming 
direction. The oblique and parallel modes are rather elec¬ 
trostatic in nature which implies that the the true nature 
of the longitudinal and transverse electric field (i.e., elec¬ 
trostatic or induced) is rather mixed and are due to sev¬ 
eral waves growing at t he s ame time, though at different 
rates (see the section |3.3| ). In order to quantify the 
role of different types of growing waves in the relativistic 
counter-steaming plasma in heating of electrons we sepa¬ 
rate (Helmholtz decomposition) the electric field into the 


rotational and compressive part, i.e. E = E'^^ such 

that = 0 and xE^ = 0. If the only waves growing 

in the counter-streaming plasma were the purely trans¬ 
verse Weibel modes then the longitudinal electric field 
would be purely inductive in nature and could solely be 
described by a rotational electric field. However, in case 



pe 


Figure 9. Temporal variation of the mean square (as in the 
figure of the electric field components E^, Ey^ and Ey 

{E = E'^ -\- E^, such that V • = 0, and V x = 0. Sub¬ 

script X and y indicate components along longitudinal (streaming 
direction) and transverse direction, respectively) are shown by the 
solid red, solid blue, dashed red, and dashed blue curves, respec¬ 
tively. During the filamentation stage the transverse electric field 
Ey is mostly electrostatic and the longitudinal electric field Ex is 
partially inductive and partially electrostatic. 


Electron 


Ion 




pe pe 


Figure 10. The left (right) panel shows the instantaneous rate 
of energy gained (lost) by the electrons (protons) due to the ro¬ 
tational and compressive components of the longitudinal electric 
field in solid red and solid blue, respectively. The dashed red and 
blue curves shows the same due to the rotational and compressive 
components of the transverse electric field Ey^ respectively. Nor¬ 
malization and computation of the mean shown here are same as 
in the figure 


of oblique and parallel electrostatic waves appearing in 
the plasma along with the purely transverse waves, the 
longitudinal electric field can not completely be described 
by a divergence free electric field. 

In figure we show to what extent the longitudinal and 
transverse components of the electric field is electrostatic 
and inductive in nature. It is evident from the figure that 
during the linear stage of the instability the transverse 
electric field is mostly electrostatic is nature and is due to 
excess of positive charges in the current filaments. The 
longitudinal electric field on the other hand is partially 
inductive and partially electrostatic, but more of an in¬ 
ductive in nature. In figure we show the work done 
by the rotational and compressive part of each compo¬ 
nents of the electric field. We find that the work done by 
the compressive and rotational part of the longitudinal 
electric field are comparable, suggesting that the role of 
resonant heating by the electrostatic waves is compara¬ 
ble with the inductive heating due to transverse Weibel 
modes. 


3.3. Parallel and oblique modes 
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Figure 11. Power spectrum for Bz, Ex, and Ey at 100 Upe 
is shown in the left, middle, and right panel panel, respectively. 
The power spectrum shows broadband nature of the Weibel insta¬ 
bility. The spectrum of Bz shows that the transverse waves are 
mostly magnetic in nature. The existence of parallel and oblique 
modes, which are rather electrostatic in nature, is apparant from 
the spectrum of Ex- 





Figure 12. The left (right) panel shows temporal evolution of 
the average amplitude of sinusoidal variation in the longitudinal 
electric (magnetic) field E^ {B^) with the wave vector being along 
the streaming direction (x-axis). The red, green, blue, and cyan 
colored curve show the field amplitude of modes with jk— 410, 
90, 16, and 8 cluope-, respectively. The average amplitude shown 
here are obtained after taking average of amplitude of any given 
mode over an ensemble of several slices parallel to the y-axis at 
equally spaced locations along the x-axis. 


We perform two dimensional spectral analysis of the 
electromagnetic field components for the simulation 
to show the simultaneous growth of several other waves 
in addition to the purely transverse Weibel waves. In 
figure we show the average amplitude of sinusoidal 
variation in electromagnetic fields along the streaming 
direction. The amplitude of variation in the magnetic 
field serves as a proxy for bending of the current filaments 
since the parallel modes are mostly electrostatic and do 
not contribute to the growth of the magnetic field. As 
evident from the figure the growth of electric field 
in the same length scale is faster than the growth of the 
magnetic field, hence by comparison a rather stronger 
electric field in waves with their wave vector along the 
longitudinal direction can be attributed to the growth 
of electrostatic modes along the streaming direction in 
addition to the contribution from the bending of current 
filaments. 

In figure we show the growth of magnetic and elec¬ 
tric field in purely transverse, purely longitudinal, and 
an oblique mode. The general trend is that the purely 
transverse modes grow the fastest and the rate of growth 
decreases with the angle between the wave vector and 
y-axis. The oblique modes, like the purely transverse 


Figure 13. The strength and sense of the magnetic field is shown 
in a small patch of the simulation box for two variant of simulations 
M64 when transverse size if the filaments is about 30 cjujp^ . Top 
panel : only counter-streaming ions beams (immobile or infinite 
mass electrons are to be assumed to be sitting at the initial loca¬ 
tions of ions in order to achieve charge neutrality). Bottom panel 
: simulation Mq 4 ^ in which case electrons are also participating in 
the plasma dynamics. Time such that the transverse size of the 
filaments is about 40 cjujpe- All modes grow faster in the absence 
of electrons which resonantly suppress the growth of waves. 


mode, are suppressed in the linear regime by the induc¬ 
tive response of the electrons which are accelerated by, 
and therefore move along the inductive electric fields. 
Note, however, that the purely transverse mode stops 
growing and at time of about 250 uOpe (figure) while the 
oblique mode continues to grow and in fact catches up 
with the purely transverse mode. This can be inter¬ 
preted as follows: The transverse mode, being the fastest 
growing, is the first to trap the electrons into the reso¬ 
nance, so that the number of resonant electrons actually 
increases. This includes those electrons that resonated 
with the oblique modes until trapped by the purely trans¬ 
verse mode. We can say that these electrons are in a 
non-linear resonance with the purely trapped mode. At 
this point, the growth of the purely transverse mode is 
impeded by the even larger number of non-linearly reso¬ 
nant electrons, so it stops growing.The oblique modes, in 
contrast, suffer less from the inductive effects of the elec¬ 
trons so they keep growing, until they have spoiled the 
alignment of the filaments made by the purely transverse 
mode. 

The role of electrons in suppressing the oblique modes 
can clearly be demonstrated by simulating the Weibel 
instability with electrons of infinite mass such that elec¬ 
trons do not interact with the plasma waves. In the ab¬ 
sence of lighter resonating electrons all modes grow faster 
and the Weibel instability give rise to rather stronger 
magnetic field (figure 13). In case of infinitely massive 
electrons the current maments which form due to the 
Weibel instability are rather short in length along the 
longitudinal direction and get quickly disoriented as com¬ 
pared to the simulation with electrons interacting with 
the waves. 

The large fluctuation in the amplitudes of the waves 
(figureis apparently due to interference of waves with 
same wave number but with different phases appearing in 































Figure 14. The left (right) panel shows temporal evolution of 
the amplitude of sinusoidal variation (obtained from two dimen¬ 
sional Fourier spectrum of electromagnetic fields) in the longitudi¬ 
nal electric (magnetic) field with wave vector being along 

x-axis, y-axis, and at an equal angle with respect to the x and 
y axes . Specifically, the red, green, and blue curves correspond 
to ( 27 t Ikx 1 Iky) = (0,100), (100,100), and(100, 0)c/c(;pe, respec¬ 
tively. The thick solid line segments in red and blue ind icate the 
growth rate predic ted by the linear theory (taken from |Shaisul-| 
|tanov et a/.|(|2012|) for electron temperature rrielmi times tlie ion 
temperature) tor the transverse and the oblique wave, respectively. 


the plasma at different spatial locations which grow inde¬ 
pendently. The phase difference in the purely transverse 
Weibel modes at different locations can also be blamed 
for the wiggle in the filaments which appears at the very 
beginning of the simulation and can contribute to the 
scattering of the electrons. 

The power spectrum of the electromagnetic fields (Fig¬ 


ure 11) reveals a broadband nature of the instabilities in 
the counter-streaming plasma. We have identified the 
growth of transverse and oblique Weibel-like modes as 
well as parallel electrostatic modes in our simulation. 
However, there can possibly be several othe r waves and 
instabilities in the counter-streaming plasma (Bret 2009), 
which make the current filaments unstable, and might 
become more pronounced for plasma parameters which 
are rather more realistic for astrophysical contexts, i.e. 
for higher ion to electrons mass ratio or larger Lorentz 
factor of the streaming plasma than simulated here. In 
any case, it is evident that whichever wave can bend the 
current filaments of streaming ions can very efficiently 
scatter electrons out of the filaments since the electrons 
have much less inertial mass as compared with ions. The 
bending of the filaments, which we suggest is due to a 
broadband nature of the Weibel instability (i.e. growth 
of the waves with their wave vector at some angle with 
respect to the transverse direction), is a crucial require¬ 
ment for the heating of the electrons as demonstrated 
by the one-dimensional simulation where poor scattering 
results in trapped electrons, which efficiently short out 
the inductive electric field created by streaming ions and 
do not get energized. 


3.4. One dimensional simulation: purely transverse 
Weibel mode 

In order to illustrate the role of oblique and electro¬ 
static modes parallel to the streaming direction we sim¬ 
ulated the counter-streaming plasma in one dimension 
such that the only modes transverse to the streaming 
direction, i.e. purely transverse Weibel models, are al¬ 
lowed to grow. The growth of oblique as well as longi¬ 
tudinal modes is suppressed by reducing the size of the 
box along the streaming direction (x-axis) to sub-skin 




Figure 15. One dimensional simulation: A) Left panel: Tempo¬ 
ral variation of the mean square (averaged over the physical do¬ 
main of the simulation) of the electromagnetic field components 
Ex (red), Ey (green), and Bz (blue) are normalized to the to¬ 
tal initial kinetic energy density nirrii + me )(70 — l)c^. B) Right 
Panel : a) Black curve : Temporal variation of Kinetic energy of 
electron KEe (normalized to the initial kinetic energy of electrons 
KEeO = (70 — l)meC^), b) Cyan Curve: Ratio of instantaneous 
Larmor radius to the instantaneous skin depth of electrons as in 
the figure 


depth level, hence reducing the simulation to effectively 
1.5 dimension. That is to say, though electromagnetic 
field and particle’s velocity are allowed to have all three 
components, they are allowed to vary along the trans¬ 
verse direction (y-axis) only. Here we discuss results from 
the simulation which has two dimensional box of size 0.3 
cjujpeX 1600 c/ujpe and 64 particles per cell. All other 
parameters for this simulation are same as in the case of 
two dimensional simulation Mq 4 ^. 

In the one dimensional simulation electrons first un¬ 
dergo the Weibel instability, which essentially is the 
Weibel instability of one species. The electron Weibel 
stage ends with the isotropisation of electrons. As the 
ions undergo the Webiel instability size of ion current 
filaments and strength of the magnetic field grow. At 
the stage when transverse size of the current filaments 
becomes ^ 2'Kclujpe electrons, which are magnetized, get 
confined in ion current filam ents between two adjacent 
peaks of the magnetic field (iWiersma and Achterberg 


2004 Lyubarsky and Eichler|2bo6 ). Electrons in the 


cur¬ 


rent 


filaments quiver in transverse direction and drift in 
the same sense as the ions. The countercurrent due to 
electrons efficiently arrest the growth of Weibel instabil¬ 
ity of the ions. Anisotropic ions continue to stream in 
the current filaments and the growth of magnetic field 
ceases (Figure [^. The lack of heating of electrons ob¬ 
served in the one dimensional case can be attributed to 
the lack of scattering of electrons out of the current fil¬ 
aments which is rather efficient in the two-dimensional 
case due to obliques waves growing along with the purely 
transverse mode. 


4. DISCUSSION AND CONCLUSIONS 

We have simulated the development of the Weibel 
instability in relativistically counter-streaming homoge¬ 
nous ion-electron beams, using the kinematic PIC 
method. The physical domain of the simulation is taken 
to be sufficiently large, in both parallel and perpendicu¬ 
lar to the plasma streaming direction, to ensure that the 
growth of Weibel unstable modes as well as the longitu¬ 
dinal electrostatic states are not frustrated due to finite 
size of the box. The homogeneous set-up simulated here 
ensures that any effect of large scale longitudinal inhomo- 
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geneity which may be present in the case of shock is sep¬ 
arated out. In the case of a shock transition there might 
be some additional heating due to cross-shock potential 
that may develop because of longitudinal separation of 
ions from electrons, since the lighter electrons can be rel¬ 
atively easily isotropised by the foreshock magnetic field, 
hence interrupting the electrons flow before the flow of 


and 


ions (iBalikhin et a/.||1993l iLyubarsky | |2006l [Lyubarsky 
i jichler||20U6| . 


Uur hndings from the plasma simulations concerning 
the heating of electrons can be summarized as follows. 


1. In the case of relativistic counter-streaming homo¬ 
geneous plasma beams electrons are accelerated to 
the energy comparable to the energy of ions and 
most of the heating occurs during the formation 
and bending of the filaments and little after the 
filaments disrupt and disorient. 

2. Comparison of the net work done by longitudinal 
and transverse electric field shows that the main 
force that takes energy aways from the ions and ac¬ 
celerates electrons is due to the longitudinal electric 
field. 


3. The work done by the transverse electric field is 
negligibly small as compared to the the longitudi¬ 
nal electric field and increasingly so for larger ions 
to electrons mass ratios (figure |^. 

4. Decomposition of the electric field into electrostatic 
and inductive component reveals that the longitu¬ 
dinal components is partially inductive, which is 
due to the growing current in the current filaments, 
and partially electrostatic, which is due to bending 
of the current filaments and growth of electrostatic 
waves along the longitudinal direction. The trans¬ 
verse component, on the other hand, is mainly elec¬ 
trostatic in nature and is due to the separation of 
charges. 

5. Calculation of rate of net work done on the elec¬ 
trons and ions by the electrostatic and inductive 
components of the longitudinal electric field shows 
that the work done by the both components are 
comparable, with that of the inductive component 
being slightly larger. It suggests that part of the 
acceleration of electrons is due to the fast moving 
chunks of current filaments (fermi-like heating, pre¬ 
sumably second order). 

6. We showed that background electrons efficiently 
take energy away from the waves generated by the 
unstable counter-streaming ions and hence signif¬ 
icantly alter the dynamics of the instabilities in 
counter-streaming plasma. It then becomes essen¬ 
tial to take continuous energization of the electrons 
into account in order to calculate growth of various 
wave modes, even in the linear approximations. 

7. Fourier spectrum of the electromagnetic fields sug¬ 
gest that waves of same wavenumber but with dif¬ 
ferent phases develop in the simulation and that 
there is significant phase mixing which can par¬ 
tially be responsible for the instability of filaments 
and scattering of electrons. 


Current filamentation and significant electron heating 
is observed in three dimensional simulations of counter¬ 
steaming plasma as well. We expect that the physical 
mechanisms in a more realistic three dimensional simu¬ 
lation should be the same as in the two dimensional case 
discussed here and can be verified in future simulations. 
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resources for the simulations. RK and DE acknowledge 
support from the Israel-U.S. Binational Science Founda¬ 
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